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ABSTRACT 

It is shown that a generahzed magneto-BernouUi mechanism can effectively 
generate high velocity ffows in the Solar chromosphere by transforming the 
plasma pressure energy into kinetic energy. It is found that at reasonable heights 
and for realistic plasma parameters, there is a precipitous pressure fall accompa- 
nied by a sharp amplification of the flow speed. 
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1. Introduction 

Recent observations, strongly fortified by immensely improved measuring and interpre- 
tive capabilities, have convincingly demonstrated that the solar corona is a highly dynamic 
arena replete with multiple-scale spatiotcmporal structures (Aschwanden et al. 2001a). A 
major new advance is the discovery that strong flows are found everywhere — in the sub- 
coronal (chromosphere) as well as in the coronal regions (see e.g. (Schrijver et al. 1999; 
Wincbarger, LeLuca and Golub 2001; Wilhelm 2001; Aschwanden et al. 2001a; Aschwanden 
et al. 2001b; Seaton et al. 2001; Winebarger et al. 2002) and references therein). Equally 
important is the growing belief and realization that the plasma flows may complement the 
abilities of the magnetic field in the creation of the amazing richness observed in the coronal 
structures. 

Interestingly enough, even before the observational mandate, theoretical efforts in har- 
nessing the plasma fiows to solve some of the riddles of solar physics had already begun. 
In particular the dynamics of fiow-based structure creation and heating was the subject of 
Mahajan et al. (1999, 2001); in this model the fiows provided the basic material as well as 
energy for the primary heating of the coronal loops. A systematic treatment of loop models 
that include flows was also developed by Orlando, Peres and Serio (1995a, 1995b), and by 
Mahajan et al. (1999, 2001). 

If fiows are to play an important and essential role in determining the dynamics and 
structure of the solar corona, we must immediately face the problem of finding sources and 
mechanisms for the creation of these flows. Catastrophic models of flow production in which 
the magnetic energy is suddenly converted into bulk kinetic energy (and thermal energy) 
are rather well-known; various forms of magnetic reconnection (flares, micro and nanoflares) 
schemes permeate the literature (see e.g. (Wilhelm 2001; Christopoilou, Georgakilas and 
Koutchmy 2001) for chromosphere up-fiow generations). A few other mechanism of this 
genre also exist: Uchida et al. (2001) proposed that the major part of the supply of energy 
and mass to the active regions of the corona may come from a dynamical leakage of magnetic 
twists produced in the subphotospheric convection layer; Ohsaki, Shatashvili, Yoshida and 
Mahajan (2001,2002) have shown how a slowly evolving closed structure (modelled as a 
double-Beltrami two-fiuid equilibrium) may experience, under appropriate conditions, a 
sudden loss of equilibrium with the initial magnetic energy appearing as the mass fiow 
energy. Another mechanism, based on loop interactions and fragmentations and explaining 
the formation of loop threads, was given in Sakai and Furusawa (2002). These mechanisms, 
though extremely interesting, are an unlikely source to account for the observed ubiquity of 
plasma fiows. One should, perhaps, look for relatively gentler, more widespread, and steadier 
mechanisms. 
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Before exploring new avenues for flow generation, and deciding which region of the 
solar atmosphere should be the first target for investigation, we seek some guidance from 
phenomenology. Based on estimates of energy fluxes required to heat the chromosphere and 
the corona, Goodman (2001) has shown that the mechanism which transports mechanical 
energy from the convection zone to the chromosphere (to sustain its heating rate) could also 
supply the energy needed to heat the corona, and accelerate the solar wind (SW). The coronal 
heating problem, in this context, is shifted to the problem of the dynamic energization of the 
chromosphere. In the latter process the role of flows is found to be critical as warranted by the 
following observations made in soft X-rays and extreme ultraviolet (EUV) wavelengths, and 
recent findings from the Transition Region and Coronal Explorer (TRACE): the overdensity 
of coronal loops, the chromospheric up-flows of heated plasma, and the localization of the 
heating function in the lower corona (e.g. Schrijver, et al. 1999; Aschwanden et al. 2001a; 
Aschwanden 2001b) and references therein). The main message then, is that the coronal 
heating problem may only be solved by including processes (including the flow dynamics) 
in the chromosphere and the transition region. The challenge, therefore, is to develop a 
semi-steady state theory of flow generation in the chromosphere. 

There are only two obvious energy sources that could power flow generation in the 
chromosphere: the magnetic field, and the thermal pressure of the plasma. We have already 
mentioned a few examples of the magnetically driven transient but sudden flow-generation. 
Looking for a more quiescent pathway, we shall now concentrate on the latter source. To 
convert thermal energy into kinetic energy, some variant of the Bernoulli mechanism (the 
existence of larger kinetic energy in regions of lower pressure) must be invoked. We will soon 
show that the double-Beltrami-Bcrnoulli states accessible to a two-fluid system in which 
the velocity field is formally treated at par with the magnetic field (Mahajan and Yoshida 
1998), can readily provide the necessary framework for this conversion. 



2. Model 

In astrophysics (particularly in the physics of the solar atmosphere) , it is useful to assign 
at least two distinct connotations to the plasma "flow": 1) the flow is a primary object 
whose dynamics bears critically on the phenomena under investigation. The problems of 
the formation and the original heating of the coronal structure, the creation of channels for 
particle escape, for instance, fall in this category, 2) The fiow is a secondary feature of the 
system, possibly created as a by-product (e.g. see Ohsaki et al. (2002)) and/or used to drive 
or suppress an instability. Since the generation of flows which will eventually create the 
coronal loops is the theme of this effort, the flows here are fundamental. 
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As pointed out earlier, our theoretical model is based on a rather simple application 
of the magnetofiuid theory developed over the last few years. We plan to restrict ourselves 
to almost steady state considerations (for a steady and continuous supply of plasma flows 
emerging from the chromosphere). Very near the photospheric surface, the influence of 
neutrals and ionization processes (and processes of flux emergence etc.) would not permit 
a quasi-equilibrium approach. A little farther distance (Ar) from the surface, however, we 
expect that the quasi-equilibrium two-fluid model will capture the essential physics of flow 
generation. 

/,From recent observational data (sec e.g. Goodman (2000); Aschwanden et al. (2001a); 
Socas-Navarro and Almeida (2002) and references therein) we could obtain the following av- 
erage plasma density and temperature at Ar ~ (500 — 2000) km: n ~ (10^^ — 10^^) cm~^; T ~ 
1 eV (for simplicity we will assume equal electron and ion temperatures) . The information 
about the magnetic fleld is a little harder to extract because of the low sensitivity and lack 
of high spatial resolution of the measurements coupled with the inhomogeneity and co- 
existence of small- and large-scale structures with different temperatures in nearby regions. 
At these distances we have different values for the network and for the internetwork fields. 
The network plasmas have typically short-scale fields in the range Bq ~ (700— 1500) G, have 
more or less uniform density and will be prone to explosive/eruptive analysis of the kind 
carried out in Ohsaki et al. (2002). The internetwork fields, on the other hand, are generally 
smaller (with some exceptions (Socas-Navarro and Almeida (2002)) — -Bo ~ 500 G, and are 
embedded in larger-scale plasma structures with inhomogeneous densities. The theory of 
steady creation of flows in the lower chromosphere will be based on these latter objects. 



To illustrate the basic physics of flow creation, we deal with the simplest two-fluid equi- 
libria with T — const leading to rT^Vp — > TV In n. The generahzation to a homentropic 
fluid with p — const • n''' is straightforward and is done in the numerical work. 

The dimensionless equations describing the model equilibrium can be read off from 
Mahajan et al. (2001) 



3. Generation of flows in lower chromosphere 




/9o In n 



) 



+ V X (V X V) = 0, 



(1) 




(2) 



V • (nV) = 0, 



(3) 
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V • b = 0, 



(4) 



where the notation is standard with the following normalizations: the density n to uq preva- 
lent at about 500 km (and farther) from the solar surface (this is where the plasma is created), 
the magnetic field to the ambient field strength at the same distance, and velocities to the 
Alfven velocity Vao- The parameters tao = GMq/V\qRq — 2/?o/tcO! «o = Ko/Re^ Po = 
c^q/V|q are defined with no, Tq, Bq. Here c^o is a sound speed, Rq is the solar radius, 
Tc — GMq/2c'^qRq, and Ajo = c/cuio is the coUisionless skin depth. 

The double-Beltrami solutions, expressing the simple physics that the electrons follow 
the field lines, while the ions, due to their inertia, follow the field lines modified by the fluid 
vorticity, are contained in the pair. 



b + CKoV X V = d n V, 



a n 



V-^Vxbl, 
n J 



(5) 



where a and d are dimensionless constants related to the two ideal invariants, the magnetic 
and the generahzed helicities (Mahajan and Yoshida 1998; Mahajan et al. 2001), 



hi 



[A ■ b) d^x, 



y (A + V) • (b + V X Y)d^x. 
On substituting (6)-(7) into (l)-(2), one obtains the BernouUi condition 

^ / 2/5orco 



/3o In n - — ) = 0, 



(6) 



(7) 



(8) 



relating the density with the flow kinetic energy, and solar gravity. Equations (1), (5), and 
(8) represent a close system, and may be easily manipulated to yield {g{r) = rca/r) 



an 



n 



-VxVxV + aoVx 



1 

an 



d nV 



+ 1 



V = 0, 



CKoV X ( -V X b ) + ao V X 



d 

an 



n — exp 



1^2 
•^0 



2g + 



+ I 1 



- b = 0, 

a 



(9) 
(10) 
(11) 



2/5o 2/3o. 

a set ready to be solved for the density, the velocity and the magnetic field. We must point out 
that this time-independent set is not suitable for studying chromospheric heating processes 
(primary heating at lower heights). The main thrust of this paper is to uncover mechanisms 
which create fiows in the chromosphere — fiows that will supply matter and energy needed to 
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create the coronal structures, and provide their primary heating. The creation and heating 
problem, of course, requires a fully time-dependent treatment (Mahajan et al. 2001). 

We have to resort to numerical methods to obtain detailed solutions for the coupled 
nonlinear system (9), (10), and (11). We have carried out a ID simulation (the relevant 
dimension being the height "Z" from the center of the Sun; Zq — Rq + Ar is the surface 
at which the boundary conditions are applied) for a variety of boundary conditions. The 
boundary surface is so chosen that at this height Zq the ionization is expected to be complete. 
For estimates of Ar given earlier, the relevant heights lie within [(1 + 0.7 • 10""^) — (1 + 2.8 • 
10-')] Rq. 

The simulation results are presented in Figs. 1-2. These are the plots of various physical 
quantities as functions of the height. The first figure consists of three frames (a-b, c-d, and 
e-f) each consisting of two pictures — one for the density and the magnetic field and the 
other for the velocity field. The parameters defining different frames are: no ~ lO^'^cm"^, 
Bo ~ 500 G, Vao ~ llOkm/s implying /?o 0.01 < 1 and = 900 for a-b; uq ~ 10^*5cm-^ 
Bq ~ 1000 G, Vao ~ 22km/s implying /5o ~ 0.2 < 1 and Vco — 950 for c-d; and Uq ~ 
lO^^m-^ and Bq ~ 1500 G, Vao ~ Hkm/s, /3o ~ 1 and = 1000 for the frame e-f. In 
each frame there are three sets of curves labelled by ckq (1-2-3 corresponding respectively to 
ao — 10~^, 10~^, 10~^), the measure of the strength of the two-fluid Hall currents. 

For all our runs the temperature at the boundary was taken to be ~ 1 eV {cgo ~ 
lOkm/s), and the boundary conditions, |6o| = 1, Vq — Ikm/s (with V^o — V^o — ^o) 
were imposed. Notice that we begin with just a small residual flow speed. The choice, 
d ~ a ~ 100 and {a — d)/tt^ ~ 10^*^ for the parameters characterizing the double Beltrami 
state, reflects the physical constraint that we are dealing with a sub-Alfvenic flow with a 
very small ao (Mahajan et al. 1999). We must admit that the values of chosen for the 
simulation are much larger than their actual values (~ 10^^ for corona and ~ 10^^^ for 
chromosphere); our present code cannot resolve the equivalent short lengths, though, we 
hope to do better in future. We believe, however, that the nature of the flnal results is 
properly captured by these artiflcial values of ckq. 

The most remarkable result of the simulation is that for small and realistic values of ao 
(curves labelled 1), there exists some height where the density begins to drop precipitously 
with a corresponding sharp rise in the flow speed. The effect is even stronger for the low beta 
(a-b are the lowest beta frames) plasmas. It is also obvious that at very short distances, the 
stratiflcation is practically due to gravity, but as we approach the velocity "blow-up" height, 
the self-consistent magneto-Bernoulli processes take over and control the density (and hence 
the velocity) stratiflcation. 
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An examination of the Bernoulli condition (11) readily yields an indirect estimate for 
the height at which the observed shock formation may take place. For a low-beta plasma, 
the sharp fall in density is expected to occur where the local flow kinetic energy exceeds the 
kinetic energy specified at the boundary (this is true for all cto), i.e, 

|V|2-T/o'>2/3o. (12) 

For the current simulation, at Po — 0.01, it occurs approximately at |V|^ > 0.02 or at 
\V\ ~ 0.14. This analytically-predicted value is very close to the simulation result (see 
Fig. 2(b)). Simulation results also confirm that the velocity blow-up distance depends mainly 
on /?o, and that the final velocity is greater for greater (3o (Fig. 2(a)). The data presented in 
Fig. 1 and Fig. 2 corresponds to a uniform temperature plasma. For this case, the variations 
in plasma pressure are entirely due to the variations in density. Since the magnetic energy 
remains practically uniform over the distance, sharp decrease in density with a corresponding 
sharp rise in the flow-speed (flow energy changes are of the order of n~^/^) is nothing but 
the expression of the commonly understood Bernoulli effect. We must emphasize that the 
general results remain unchanged in our extensive simulations in which the temperature is 
allowed to vary (but we have to use a homentropic equation of state to analytically derive 
the Beltrami states. The final parameters, naturally, depend upon the adiabaticity index 7). 

Taking into account the fact that the fiducial height Zq is different for different cases 
(larger for smaller /3o), one expects that the plots for the "blow-up" distance will approx- 
imately match one other, and there exists some minimum distance (possibly of the order 
500 km or so) from the photosphere below which steady up-fiows can not be generated. This 
expectation is simply supported by the numerical results given in Fig. 2(a). Observations 
clearly indicate that there exists a narrow layer above the photosphere where no spicules, 
mottles etc. are seen. 

To check whether the generated flows arc predominantly radial or somewhat more 
isotropic (to explain the observational constraints) we studied in detail the relatively large 
/3o case (fixing j3o is quite difficult due to complications like ionization) and found that the 
fiows tend to be mostly radial only for large ao (see, for example, plots labeled 2 and 3 in 
Fig. l(b,d,f)). The situation could change considerably when we deal with a more inclusive 
time-dependent dynamical model with dissipation. Plasma heating, then, could result from 
the dissipation of the perpendicular energy so that at larger distances, the flows would have 
larger radial components. Heating would also keep /3(r, t) large at upper heights shifting the 
velocity blow-up distance further or eliminating it all together; we know from Fig. 1 that as 
/9o goes up, the density fall(vclocity ampliflcation) becomes smoother. These issues will be 
dealt with later in a more detailed work. 

Note, that if one were to ignore the flow term in (8) (a totally wrong assumption com- 
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monly used in many studies), we will end up finding essentially radial flows. The magnitude 
of these flows, however, remains small; there is no region of sharp rise (12), and the generated 
flows achieve reasonable energies at heights typically 10 times greater than the heights at 
which the correct Bernoulli condition would do the trick. 



4. Conclusions and Summary 

We have shown a possible pathway for a steady generation of flows in the quasi- 
equilibrium stage established through the interaction of the emerging magnetic fluxes with 
the existing cold solar chromosphere (when an ionized ~ 1 eV plasma is trapped in in- 
ternetwork structures). The suggested mechanism is a straightforward application of the 
recently-developed magnctofluid model (Mahajan and Yoshida 1998; Mahajan et al. 1999, 
2001); a generalized Bernoulli mechanism (a necessary condition for the double-Beltrami 
magnetofluid equilibrium) allows the pressure energy to be very effectively transformed to 
flow kinetic energy as the plasma moves away from the sun. We flnd that at reasonable 
heights and for realistic plasma parameters, there is a precipitous pressure fall with a sharp 
amplification of the fiow speed. In the presence of dissipation, these flows are likely to play a 
fundamental role in the heating of the inner and upper chromosphere, although our explicit 
purpose in this paper was to create a steady source of matter and energy for the forma- 
tion and primary heating of the corona. Our preliminary results agree with the observation 
data, and lend promise to attempts, based on the existence of subcoronal flows, to tackle 
unresolved problems like the coronal heating and origin of the solar wind. 
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Fig. 1. — Plots of density, magnetic fields and velocity versus height for values of ao and f3o. 
Sub-figures (a) and (b) are for Pq = 0.01, = 900; (c) and (d) are for /3o = 0.2, = 950; 
(e) and (f) are for Pq = 1, = 1000. The numbers 1, 2, 3 represent ao = 10~^, 10~^ 10~^ 
respectively. Vy is not displayed since its behavior is practically similar to Vx- The velocity 
blow-up is controlled by f3o. For a bigger (unrealistic) ao there is a splitting of the velocity 
components — at the end the radial component is dominant. Magnetic field energy does not 
change much on these distances. 

Fig. 2. — The "blow-up" distance (a) and velocity (b) versus ao. The smaller the f3o, the 
smaller the "blow-up" distance and smaller the velocity at "blow-up" (compare with (12)). 
For fixed Pq, the process is less sensitive to changes in aQ. 



